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Abstract 

We calculate the far-field stress-energy tensor of a quark gluon plasma due to a heavy- 
quark meson moving through it, using a semi-classical description in string theory of the 
meson as a string hanging into anti-de Sitter space. We find that these mesons create a 
shock wave but no diffusion wake, in contrast to single heavy quarks described via a trailing 
string, where the diffusion wake is strong. If thermal QCD responds similarly to heavy 
quarks and heavy quarkonia, the presence or suppression of the diffusion wake constitutes a 
prediction which can in principle be checked in heavy-ion collisions. 
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1 Introduction 

As a quark moves through a quark gluon plasma (QGP), it loses energy to the medium. An 
extensive theoretical enterprise aims at describing this energy loss using some combination 
of perturbative QCD methods and hydrodynamics: for recent reviews see for example [HIS]. 
There is also a substantial effort to compare energy loss in QCD to calculations in gauge 
theories with known holographic duals, especially A/" = 4 super- Yang-Mills: papers in this 
genre that directly set the stage for the current work are [3], HI O [6], [3, [H [9], [TOl [11], |12], and 
related works include [131 [IH [13 [IS]- In [12] it was shown that in a large class of theories, 
whose gravity duals may be described by asymptotically AdS black hole geometries, the ratio 
of energy from an infinitely heavy quark going into sound modes to the energy coming in 
from the diffusion wake behind the moving quark is 1 + f ^ : 1. Thus the diffusion wake seems 
to be an important and persistent feature of the dynamics. On the other hand, according to 
the analyses of [I71[T8], PHENIX data supporting the hypothesis of jet splitting [19] appears 
to favor a scenario where the wake is either absent or relatively suppressed. It is worth noting 
that published results from STAR [20] with more inclusive transverse momentum cuts do 
not show jet splitting so much as substantial medium-induced jet broadening. 

In this context it is interesting to inquire what pattern of energy fiow surrounds other 
probes of A/" = 4 super- Yang-Mills theory. Next to a single infinitely massive quark, described 
by a string with one end attached to the boundary of AdS-Schwarzschild and the other 
trailing behind it, the simplest such object is a heavy quarkonium-like state, described by a 
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string with both ends attached to the boundary [21]. Such objects have been considered in 
[H [H] , where it was observed that they experience no drag force, provided their velocity is 
not too large relative to their size times the ambient temperature. One may speculate that 
the absence of a diffusion wake is generally associated with no-drag configurations, meaning 
configurations where the drag force is not visible in the approximation used to describe it 
(large N in our case)ll] 

The aim of this paper is to compute the large distance asymptotics of the gauge theory 
stress-energy tensor produced by such heavy quarkonium states, both when the velocity is 
parallel to the separation of the quark and anti-quark, and when the velocity is perpendicular 
to it. We find that in both cases a sonic boom is observed but no diffusion wake. This may 
be compared with a phenomenological model [22] used to explain the jet-splitting effect. 
There it was assumed that when a single quark moves through the medium, the total energy 
loss is a subleading effect and the diffusion wake is suppressed. Since [TTl [T2] show that this 
assumption does not hold for heavy quarks in a large class of theories with asymptotically 
AdS duals, our result provides a first hint that jet-splitting (to the extent that it really 
occurs) may involve the motion of composite or colorless particles through the QGP. 

There are two main steps in the computation of the gauge theory stress-energy ten- 
sor. First, in section [21 we review the string configuration of interest and identify its 
five-dimensional stress-energy tensor. Second, in section [3|, we solve the linearized five- 
dimensional Einstein equations sourced by this five-dimensional stress-energy tensor. Our 
solution is entirely analytic but is only accurate to leading order in small wave-numbers. 
From it, the long-distance behavior of the gauge theory stress-energy tensor can be read off 
directly. We continue in section |4] with some comparisons to hydrodynamics and a discussion 
of the maximum velocity for the mesons we are considering. 



2 Test string description of heavy-quark mesons 

The purpose of this section is to review the string configurations that describe the heavy- 
quark mesons of interest and to find out how these configurations source perturbations of the 
bulk metric. The starting point is the Einstein-Hilbert plus Nambu-Goto action, conveniently 
expressed as 

^S.S.G. thanks H. Liu for a discussion that suggested this speculation. 
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Here G^j, is the space time metric, and gab is the induced metric on the worldsheet. The X'^ 
are the embedding functions for the string. The background hne element is 
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h{z)dt' + Y^{dx'f + 



Kz) 



(2) 



i=l 



where 



L 



z 



A 



and 



h{z) = 1 




(3) 



z 



This translation-invariant black hole background describes an infinite, static, thermal medium 
of A/" = 4 super- Yang-Mills. 

For simplicity, we will only consider the cases where the quark-anti-quark separation is 
parallel to the direction of motion, or perpendicular to it. In a gauge where = z and 
X° = t we may use the ansatze 



to describe the shape of the string [3 [H [HI [15] . We shall usually drop the || and ± subscripts 
below, re-introducing them only when necessary. One should be able to generalize the results 
here to the less symmetric case where the meson separation is at an arbitrary angle Q relative 
to its velocity. It would also be interesting to consider extending this analysis to finite quark 
masses and to rotating mesons, along the lines of [H]. 

With the ansatz is double valued, so we need to specify ^ separately in the range 

— vt>0 and X^ — vt<0 for the parallel case, and in the range X"^ > and X^ < in the 
perpendicular case. Consider first the solution at X^ — vt > (or X^ > 0). In this region, 
we define ^+ = ^ — i being the distance between the quarks. If z^ is the maximum value 
of z along the string, then the boundary conditions satisfied by ^+ are 



Similarly, at X^ — vt < (or X"^ < 0) we define = C, + |£, which satisfies the boundary 



= {t,vt + ^ii{z),0,0,z) 
X^ = {t,vt,^^{z),0,z) 



(4a) 
(4b) 




(5) 



conditions 



e-(o) = o. 



(6) 
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Putting these two solutions together, one finds a continuous, different iable curve describing 
the shape of the string. 

Because the Lagrangian is time independent, the momentum 11^ conjugate to ^ is con- 
served. Exphcitly, these momenta are 



1 CM^ 

n.„ = , " (7a) 



We define the quantities 

j-j ^ V^(^m)a(^m)^ X 

= W ' ^^^^ 

which coincide with (I7j) for ^ = ^+ after implementing the boundary condition ^^(2;m) = — oo. 
By inverting ([7]), ^+ can be written in integral form. For the case where the string is 
perpendicular to the direction of the velocity, this integral can be carried out explicitly in 
terms of Appell functions [23]. From symmetry, we have = Also i = —2^^{zm)- Of 
the two possible configurations which solve i = —2^+{zm) only the one in which the string 
is closer to the boundary is stable [231 1211 125]. We shall work with this configuration. 

By definition, the bulk stress-energy tensor tensor of the string is the quantity J^j^ that 
enters into the Einstein equations as 



1 6 2 



Starting from ([T]), one straightforwardly finds 



= —^^dX^^dX^Six^ -vt- X^)5{x^ - X^)5{x^)Q{zm - z) (10) 

2Txa' V— G 



where 



,1 a; > 

e(x) = { (11) 
a; < 0. 

This stress-energy tensor may conveniently be expressed in terms of a co-moving Fourier 
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transform, meaning a three-dimensional Fourier transform in the x directions, with time 
dependence included by using — vt in place of x^ in the exponent: 



(2^ 



Denoting 



we find 



(^) 



a = cos(fc, (e+ + ^/2)) 
5, = sin(A;, ({+ + £/2)) , 



n,. 



2Civ/hi'^ 



y 2iSiv 
for the parallel configuration and 



n, 



-Q{Zm - Z)X 



( 2C2(1 + K'^)/hi'_, 2C2V{1 + hi'^)/hi'^ 

2C2v{l + hC^)/he+ 2C2v\l + H'^)lh^'^ 





y 








(12) 



(13) 







2Civ/h^'^ 
-2Ci{-v^ + h^^'^)/h^'^ 



2iSih 



2iSiv \ 
2iSih 



2Ci{v^ -h)/i'^) 








(14) 













2C2{v^-h)^'^ -2iS2{v^-h) 



-2iS2{v^-h) 2C2{v^-h)l^'^l 



(15) 



for the perpendicular one. 



3 The gauge theory stress-energy tensor 

The computation of the gauge theory stress-energy tensor proceeds by solving the linearized 
Einstein equations sourced by fll4p or (1151) . This computation is non-trivial even in the 
leading long-distance approximation we will use in section 13.31 To set the stage, we review 
in section 13.11 the general approach for computing the holographic stress-energy tensor, and 
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we discuss in section [3^ the symmetries of the problem as well as the counting of integration 
constants. 

3.1 Computing the holographic stress-energy tensor 

The metric perturbations sourced by the mesonic string are dual to fluctuations of the stress- 
energy tensor of the plasma according to the usual dictionary of AdS/CFT [26l [27] . Further 
developments, for example [28l [29l [30] , have clarified how to use this dictionary to perform 
calculations of appropriately regularized one-point functions of the type we are interested in. 
The approach used here is essentially that of [6], which gives the same results in this case as 
the more formal procedures of [291 EQ] once certain divergent contact terms, to be described 
below, are subtracted. We start with a metric of the form 

G^^u = Gjl] + «:5«'^/.. , (16) 

where G^^} is the AdSs-Schwarzschild background ([2]) and H^,, is small compared to 

(at least close to the conformal boundary.) The procedure is first to pass to axial gauge, 

Hp,z = 0, and then to expand the remaining components near the conformal boundary: 

where we have already used the boundary condition HmniO) = implying that the boundary 
theory metric is Minkowski. Without this boundary condition we would have gotten 

Hmn ~ Rmn + SmnZ + (D{z ). (18) 

The indices m and n run from through 3 corresponding to the Minkowski space coordinates 
{t,x^,x'^,x^). With the boundary conditions leading to (fTTj) . the leading behavior of Hmn is 
PmnZ^, but the coefficients Pmn are found to be analytic in k, corresponding to singularities 
with delta-function support in position space. In principle, these singularities should con- 
tribute to the stress-energy tensor, and their contribution includes an infinite prefactor due 
to the scaling with rather than z^. Such divergent terms also appear when considering 
static quarks, and they are a result of the divergent mass of the quarks. In |2T] it was 
argued that these types of contributions should be dropped in order to obtain finite energy 
configurations. As in [6], we shall use a prescription where these terms are subtracted. A 
more detailed analysis along the lines of [291 EO] should give identical results. With this 
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prescription, the expectation value of the gauge theory stress-energy tensor is 

^y^mn) (-^m.n) bath + (5T„„) , (19) 

where 

(T„„)bath = ^iV'T^ diag{3, 1,1,1} (20) 

o 

and 

(5T„„) = 2L'Qmn (21) 

is the contribution from the meson. The coefficients Qmn are non-analytic in fc, so they 
capture long-distance behavior of {6Tmn)- As we will see, they signal the existence of a sonic 
boom but no diffusion wake. 



3.2 Symmetries and the counting of integration constants 

In principle, the analysis ending in fl2T]) can be used to extract the holographic stress-energy 
tensor even for metrics that are finite deformations of the AdSs-Schwarzschild metric: all 
that matters is that the deformation should become small near the conformal boundary. 
But it is very difficult to solve the full non-linear Einstein equations in the presence of a 
localized source like the strings described in section [2l We therefore pass to the linearized 
approximation of the Einstein equations: 

V^^H,^ = J^, , (22) 

where 1^^^^^ is a variant of the Lichnerowicz operator. The linearized approximation should 
be valid at large N, because in this approximation, the radius of AdSs is much larger than 
the length scale in five dimensions where gravitational fields sourced by the string become 
strong. 

Recall that we work throughout in axial gauge, where H^^ = 0. Also, we use the co- 
moving ansatz for metric perturbations, 

HUt, x,z) = J ^e^[^i(-^-*)+^^-^+'=3x3]^_(^^ ^) . (23) 

Plugging ([23]) into 0221) . with either (IT^ or fllSp as source terms, yields fifteen inhomogeneous 
differential equations for the co-moving Fourier coefficients Hmn{k,z). These are ordinary 
differential equations in the variable z. An extended discussion of their structure and the 
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appropriate boundary conditions in the case of a cylindrically symmetric source (such as the 
inhne meson) was given in [B]. 

The ten equations with /i and v ranging from to 3 (briefly, the mn equations) are second 
order, so it takes twenty integration constants to specify a solution. These twenty constants 
are the coefficients Rmn and Qmn entering into f|T7|) and f|T8l) PI and setting Rmn = amounts 
to fixing ten of them. The five equations with fi = z are first order constraints that fix five 
additional integration constants. The constraint from the zz equation is 

Qoo - Qu - Q22 - Q33 = , (24) 

which is nothing but the condition for tracelessness of the stress-energy tensor of the A/" = 4 
theory. If we denote 

k"" = {vk^ k, k2 fcs) , (25) 
then the four remaining constraints, from the mz equations, are 

ik^Q^^ = iL-^Ii\\sxYi(^^{v -100) (26a) 
for the inline configuration, and 

iA;"^Qnm = sin (^^^ (0 -1 o) (26b) 

for the perpendicular configuration. Equations fl26p represent the energy conservation equa- 
tions in the boundary gauge theory. The remaining five integration constants are determined 
by the boundary conditions at the horizon which we discuss below. 

3.3 Solving the linearized Einstein equations for small k 

Having summarized the information that can be extracted from the first order constraint 
equations (the nz equation in the axial gauge notation of the previous section), it is now time 
to attack the second order equations, coming from mn components of the linearized Einstein 
equations fl22|) . It was shown in [6] how cylindrical symmetry around the direction of motion 
allows these equations to be partially decoupled; however, this method helps us directly only 
in the case where the quark separation is along the direction of motion. The key simplification 

^The Pmn are uniquely fixed by the behavior of the string near the boundary. 
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arises from isolating the leading small k behaviorjfl This behavior is expected to capture the 
large x asymptotic behavior of {6Tmn) because the propagating modes are hydrodynamic: 
non-hydrodynamical modes have some attenuation length and therefore decay exponentially 
at large distances from the source. The simplification at small k is that the mn equations 
naturally decouple and take the form 

""'^""t- (^'«'^") + Oi^) = Sx , (27) 



where X corresponds to various linear combinations of Hmn and Sx is the corresponding 
combination of Jmn'sE More explicitly, we find that Hqi, Hq2 and i/os satisfy (1271) with 
n = and 

2 

Soi = —-^Jqi ■ (28) 
The other off-diagonal components i/12, i^is and /f23 satisfy ( l27l) with n = 1 and 

2 

Sij = ■ (29) 

The diagonal components need to be reshuffled in order to be brought in the form ( 127|) . 
Defining 

-fo = Hii + H22 + H^3 

Fi = --H11 + -H22 + -i^33 

i 6 i (30) 

F2 = —H22 + -f^as 

-^3 = + Hii + H22 + -f^33 

h 

we find that each of the Fi satisfies equation (1271) with no = 1/2, rii = 77,2 = 1 and 773 = 3/2, 
and ^ 

= ~j^Joo 

Si = (-2Jn + J22 + J33) 

2 (31) 
'S'2 — — (^33 — J22) 

_2 fl 

'S'3 — — I ■^■^00 + ^11 + J22 + J33 



•^Short distance (large k) asymptoties are likely to be obtained via the methods developed in [311 [32l [33] . 
"'With a slight abuse of notation, we shall only include the indices of X in the subscript of Sx- Thus, in 
([28]) and ([29]) . instead of Sh^„ we shall write Smn- Similarly, in ([50)) and ([?T|) we use Si instead of Spi ■ 
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The solution to equation (!27|) is 



a{zfh{zrX\z) = qx + r'- K ) K ) xy ) m ^g^) 

where the two branches of the solution arise from the step function in ( fTOj) . Thus, 

=r,+ dy^4^ +r ' -'^''''^'^ (33) 

Jo a{yYh{yY ^ ^ , 

where gx and are integration constants. We will find it convenient to define 

T.x{z)= I dxh{xYa{xfSx{x). (34) 

•J Zj-ir 



Comparing fl33l) to flTTl) . it is clear that the are related to Rmn, and the are related 
to the Qmn- So we should set rx = for all ten choices of X in order not to introduce a 
deformation of the metric in the boundary theory. 

Even before working out the precise relationship between the qx and the Qmn (which 
follows from asymptotic expansions in small z and use of (l30l) ). we reason that the constraints 
(!24|) and (126|) must translate into five constraints among the qx- In order to uniquely specify 
the physical solution to the linearized Einstein equations, five more boundary conditions are 
needed. As in the case of the trailing string, these come from horizon boundary conditions, 
which demand that there should be no outgoing modes at 2; = ^h- Let's focus (slightly 
presciently) on the five choices of X where n = 1. Then we see from fl33l) that, near z = z-n, 

X{z) = -——In h{z) + dy 

Jo (^{yfKy) (35) 

~ qx \og{z^ - z) + ... . 

The claim is that these qx (i.e. the five corresponding to the choices of X where n = 1) 
have to be zero. To see this, we reason that at finite /c, outgoing modes have the form 
X{z) ~ (^H — 2;)™'''i^h/4 j^gg^j, 2; = z-a., while infalling ones behave as X{z) ~ (^h — 2;)~*''^i^h/4_ 
At finite /c, we would allow the infalling solution and disallow the outgoing one. In the limit 
of small /c, one may expand 

(.H - ^)±-'=-H/4 ^ 1 ± l^g(^^ _ ^) _ (36) 
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The coefficient of the log term in (136|) is suppressed by a power of k, hence the claim that 
gx = in What is perhaps surprising about the argument following (15^ is that, at 

leading order in k, it doesn't distinguish between infalling and outgoing modes. But that 
is in fact appropriate: dissipative effects, such as shear viscosity, are invisible at this order 
in k, so it's right to demand that there is neither inward nor outward flux at the horizon. 
When a similar analysis is carried out for the three choices of X where n = and the ones 
with n = 1/2 and n = 3/2, the outcome is that there is no constraint from the horizon on 
the corresponding qx- 

Having argued that all ten rx and five of the qx must equal zero, then together with the 
five constraints fl24l) and (1261) we can determine all the Qmn- For the choices of X where 
n = 1 we write 

X{z) = Pxz^ + Qxz'' + 0{z''), (37) 

where Px is related to the Pmn terms in (fT7|l and Qx to Qmn- Since T,x{z) = 0{z~^), it 
follows from (1331) that 

Qx = T^lim9,(^Sx(^)) . (38) 

For both the inline and perpendicular configurations, the above equation combined with (l34l) 
and the corresponding source terms of the Qu, Q23, and Qis equations give 

Q12 = Q23 = Qi3 = 0{e) . (39) 



Thus far, the analysis is essentially the same for the perpendicular and inline cases. In the 
remainder of the analysis we treat the perpendicular and inline cases separately. 

For the perpendicular case, flMl) combined with the appropriate source terms of the Fi 
and F2 equations gives 

si(.) = --n^ Liz) + m - r dx ) 

3 V Jzm [v^-h)hi^J (40) 

T.2{z)=m^ (e+(^)+£/2) 

where we have used ^+(2;^) = — ^/2. The integral in (HOj) and ^+{z) may be written in terms 

finer analysis shows tliat 513 = (723 = 912 = 0{k^) while qi — 0{k) and q2 — 0{k'^) for the inline 
meson and (72 = 0{k) for the perpendicular one. 



11 



of Appell hypergeometric functions. From (138|) and (HOj) we find 

1 



Qi 



6L3 



n, 



(41) 



2L3 



n, 



where 



91/2 3/2 r 1 1 4i 

Z TT Z„ 2-ri 4' 2' 4' 



(42) 



At this point, we have determined five of the Qmn for the perpendicular case. The rest follow 
from conservation and tracelessness of the stress-energy tensor (or, in gravitational terms, 
the first order constraints). Explicitly, by inverting the relations (130|) and using (!24l) . (I26bp . 
( !39|) and (14T|) we can obtain the stress-energy tensor through ( !2T|) . and it reads as follows: 



n 



(43) 



with 



(^) 



f —k"^ k\v kik2V kiksv^ 

kfv -klv^ 

kik2V 2v'^kl - A;2 q 

\A:iA;3t; -kjv^J 



(44) 



and 



/ t;(3A;2-A;2) 
fci - 2ky 
-k^k2 (1 - 



ki - 2ky 
-V {ki - 2ky) 





-k,k2 (1 - V') 



klv{l-v^) 




-A;iA;3 (1 - i;2)\ 



klv{l-v^) J 



(45) 



For the inline case, we find that S2 = 0{k) so that Q2 = 0{k), but 



(46) 



Here we cannot carry out the second integral explicitly. From fl38|) it follows that 

n,. 



Qi 



3L3 



(£ + 2t;V||) 



(47) 
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with 



z-^o dz 



dx ., , , (48) 
which can be written as a definite integral and carried out numerically, if desired, for specified 



V and Zm- Thus, 



where 



n,. 



/c^ — 3klv'^ 



(^r(il)+2w||rtl,) 



(49) 



(II) 



/-3i;2A;2 - P 
2klv 
2kik2V 
\ 2kik^v 



2klv 
v^kj - e 





2kik2V 


-2klv^ 




2kik^v \ 



-2klvy 



and Tmn is the same as in fj^5l) . A mild consistency check of these results is that in the v 
limit, Tj_mn coincides with T\\mn up to a ninety-degree rotation. 

Equations P3|) - fH5|) and fH9|) - fl50|) are our final results from solving the linearized Einstein 
equations in the presence of the mesonic string configurations, at leading order in small k. 
The two poles, located along the real ki axis at k\ = k'^/{3v^ — 1) for t;^ > 1/3 indicate that 
a shock wave will appear in real space. The absence of a pole at /ci = indicates that there 
is no diffusion wake [9]. 



(50) 







4 Discussion 

Calculations in AdS/CFT at small wave- number may reasonably be expected to relate to 
hydrodynamics. This theme has been emphasized, for example, in the recent review [M] ; see 
also [35] and references therein. We pursue this comparison in section 14. 1[ In section 14.21 we 
remark briefly on the well-known speed limit for flnite-mass mesons constructed as hanging 
strings and make some tentative remarks about the possible phenomenological relevance of 
our calculations. 
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4.1 Comparison to hydrodynamics 

From (12^ it is clear that the source for the mesonic configuration has a dipole structure. In 
real space (at t = 0) equations fl26|) read 

d^T\\^,^ = {v, -1, 0, 0)„ n,! - £/2) - + 1/2)) 5{x2)6{x^) (51a) 
d'^T^mn = (0, 0, -1, 0)„ n^(5(xi) {5{x2 - ill) - 6ix2 + i/2)) 6{x-s) . (51b) 

Let us compare flMl) to a hydrodynamic analysis. In the linearized hydrodynamic approxi- 
mation, the stress-energy tensor may be written in terms of the velocity field. This amounts 
to six relations 

^hydro ~ ^^^^''T^yjj.Q — -T s (^^T^^ydio + ^^"'^hydro ~ g ^^''^'^hydro^ (52) 

between the space-space components of the stress-energy tensor, the energy density, and the 
Poynting vector. The four additional relations needed to completely specify the stress-energy 
tensor are 

if^mT^ydro — Aydro • (53) 

Naively one would like to identify ./hydro with the source terms in (I5T1) . However, in [22l [T2] 
it was argued that, in the case of a moving quark, the quark together with its near field 
source the hydrodynamic modes, as opposed to just the quark itself. 

In [12] a (non-unique) prescription was used to compare the full stress-energy tensor, in 
our case P3|l and pojl . to the one expected based on a hydrodynamic analysis: by identifying 

(5r™°) = T-°„ , (54) 

and using (l52l) we find 

(^Tii) =T||hydro + ^n|| (£ + 2i;V)diag(o -2 1 l) (55a) 
(5Tx) =T^hydro + ^nVc^diag(o -2 1 l) +^n^£diag(^0 1 -2 l) . (55b) 

With these definitions, the hydrodynamic sources are given by 

Jl^y^^^ = -U (^ikiviS liki{i-Av^a) \ik2{i + 2v^a) \ik^{i + 2v^-a)^ (56) 
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where 6 = 1 for the parallel configuration and zero for the perpendicular one. Adopting the 
notation of [22] where 

Jn = (eo go + ki9i k2gi k^g^ (57) 

we find 



e\\ = ikiUvi 
go = 2ikillv'^a 



e_Lo = 



(58a) 
(58b) 



The source (1571) was shown to produce a diffusion wake when gi vanishes, and following 
a Cooper- Frye treatment of hadronization the authors of [22] concluded that no jet-splitting 
was predicted. A similar analysis with g^ = eo = led to a prediction of jet-splitting, 
although subsequent work [361 [H] suggests that the rate of parton energy loss may have to 
be tuned to unrealistically high values in order to match to data. In |12] a figure of merit 71 
was proposed to distinguish quantitatively between the two scenarios: 



gi{k=l/R) 



eoik = 1/R) 



(59) 



The expectation is that large 71 leads to a prediction of jet-splitting, while small 71 does 
not. Since e_Lo = we find immediately that 7i_l ^ 00, though perhaps a more conservative 
approach would be to redefine 



gi{k=l/R) 



which gives 



71 



goik = 1/R) 



1 



(60) 



(61) 



Caution is in order when using (158|) . fl59l) . or fl60l) . because it's not clear that the quantities 
Co, go, and gi can be defined in a way that is independent of contact terms, meaning terms 
which are analytic in k. The reason for this is that while {6Tmn) has non-analytic behavior 
in k and can therefore be understood as encoding universal behavior about the long-distance 
asymptotics, this non-analyticity vanishes when we use (!53|) and (15^ . and some prescription 
for the treatment of contact terms is implied. Despite these intricacies, the main qualitative 
features of our analysis is easy to state: for the heavy meson there is a sonic boom and no 
diffusion wake, and it's because of a dipole structure of the source. 
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4.2 Restrictions on the velocity 

It is tantalizing to have a microscopic description — a rather simple one at that — of a source 
that avoids producing a diffusion wake. But is this source just a contrivance, or does it 
have some relevance to QGP physics? The charm quark is heavy enough compared to the 
temperature scale at RHIC for the calculations we have presented to be compared meaning- 
fully to the behavior of the J / tp meson moving in a quark-gluon plasma. The J/ ip is indeed 
produced in heavy-ion collisions at RHIC, and some of its properties are studied: see for 
example [37] as well as results from CERN experiments at lower energies, e.g. [38] • If one 
could construct a di-hadron correlator where one hadron is tagged as coming from a J/ip, 
it would be reasonable from our analysis to expect a split-jet structure, corresponding to a 
suppressed diffusion wake. But J/V' production is so rare at RHIC that it can't possibly 
be relevant to the medium modification of the away-side jet as seen in untagged di-hadron 
correlators, studied for example in [2DI [TU] . In these studies, the overwhelming majority 
of energetic partons are probably energetic light quarks and gluons — or at least, states un- 
related to heavy flavor. But "light" at RHIC energies probably should be understood to 
include the strange quark. It is tempting then to speculate that strange quark mesons, for 
example the 0, might exhibit the effect we describe, namely production of a sonic boom but 
not a diffusion wake (or a suppressed diffusion wake). Even if it seems improbable that such 
an effect could explain jet-splitting as observed in [19], one might reasonably inquire whether 
strange quark mesons are enhanced in a directional fashion in the final state. 

One immediate question that comes to mind when considering this scenario is how the 
(f) forms soon enough to have time to propagate through a significant amount of quark- 
gluon plasma. With our current techniques, our AdS/CFT analysis does shed light on this 
problem. What we do note, is that there is an upper bound on the energy of the meson 
which follows from a bound on its velocity; finite mass quarks described by strings ending 
on branes (and mesons made from such quarks) cannot move faster than a certain velocity, 
calculated from the condition that the the endpoints of the string cannot move faster than 
light. This observation seems to have been made in various nearly equivalent ways by several 
authors, for example [HI [TJ El EH]- Here we will follow approximately the discussion of [39] , 
adapted to the case of mesons. 

The setup for finite-mass quarks involves a D-brane located at a finite value ^ = ^* in 
the AdSs-Schwarzschild background, our equation ([2])|^ The thermal mass of a static quark, 

^Most simply, this can be a D7-brane, as in [40]. Note that in this case, what means is the deepest 
inside AdSs that the D7-brane reaches. 
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first calculated in [3], is 

An endpoint of the string moving at constant depth z^, with a fixed coordinate velocity, say 
= vt, has an induced line element 

^^induccd = [-h{z.) + v'] de . (63) 

Evidently, this induced line element is timelike when f < w^,, null when v = v^, and spacelike 
when V > v^, where 

v, = ./h{^. (64) 

Thus we conclude that is the maximum velocity that the meson is allowed to have. The 
associated maximum Lorentz gamma factor has a simple form: 



2 



where in the last equality we have remembered that T = 1/ tt^h and Oym-^ ~ -^^/ Without 
attempting to account for the possibility of non-trivial dispersion relations, the maximum 
energy of a meson built as a string with both ends on the D-brane at z = is 

^ =^rn =m 11+ '^^^^^^^ ] (66) 

-"-^max /*"''mcson '""meson I „ / — I • V / 

As a simple-minded example, imagine the meson, with ^ 1 GeV, to be composed 
of s and s with ~ 300 MeV. Suppose this meson propagates through a plasma with 
T = 250 MeV and g^^^ = 10. Then E^^ax ~ 3 GeV, and ^ 0.95. While is comfortably 
above the speed of sound, -Emax is slightly below the energy range in which there is some 
experimental support for the phenomenon of away-side "punch-through" (see for example 
[^)FI One should keep in mind that in this example we've been applying an A/" = 4 result 
to QCD, so that a precise bound on the meson mass cannot be obtained. It should be clear 
though, that decreasing the quark mass in (166!) brings -Emax down quickly, especially if mmoson 
is also decreased. Thus, a heavy-light meson, which would be described in this context by a 



^It should be noted that the peak at A<f> = tt for high pT particles is not unambiguous evidence for 
high-energy particles going all the way through the QGP: an alternative scenario hinging on surface bias 
may also explain the data. S. Gubser thanks P. Steinberg for a discussion on this point. 
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straight string stretching between two D7-branes [l2l [3l HS] , is also problematic. 

Our bottom line is that, while single heavy quarks as described in AdS/CFT generate a 
strong diffusion wake in addition to a sonic boom, heavy quarkonia in the same description 
generate only a sonic boom. Applying these heavy quarkonium results to the cf) meson is 
risky. Even so, it would be interesting to inquire whether there is any directed excess of 
strange mesons associated with hard probes. If such an excess exists, how does it relate to 
jet-splitting? 
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